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FOREWORD 


The purpose of this compilation of theoretical results is to provide an 
internally consistent set of model atmospheres that can be used for 
comparison with observations. Enough data on the depth-dependent quan- 
tities are given so that other quantities can be easily computed. The 
observable quantities—monochromatic emergent fluxes, discontinuities, con- 
tinuum slopes, UBV and uvby colors, and equivalent width and central 
intensity of H^-~are presented. 
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BASIC EQUATIONS 

The models were computed with standard techniques that have been 
described elsewhere (Mihalas, 1967; Pecker, 1965). The actual computer 
programs used have been described in detail by Klinglesmith (1967). 

The equation of hydrostatic equilibrium included a radiation pressure term 
(Woolley & Stibbs, 1953) and was solved with a fourth-order Runge-Kutta 
method with variable step size. The monochromatic equation of transfer in a 
plane-parallel semi-infinite atmosphere was solved with the method of 
parabolas (Avrett and Loeser, 1963), and the models were flux corrected 
through the use of Lucy’s temperature correction procedure (1964). 

The continuum opacity sources used are listed in Table 1 with the 
references from which the equations were taken. The approximate expres- 
sions derived by Griem (1960) and modified by Underhill (1962) were used 
to compute the bound-bound opacity of the Lyman and Balmer series of 
hydrogen. The total line opacity at a particular wavelength is a sum of aU the 
individual line opacities with upper-level quantum numbers less than or 
equal to 40. Near the series limits, a photon may be absorbed into either 
a continuum or line transition. This probability (X) is described by 
de Jager and Neven (1960), and Mihalas (1966) shows how it can be used 
to take into account, in an approximate manner, the influence of the 
confluence of the Balmer lines. 


TABLE OF OBSERVABLES 

Table 2 contains the observable quantities for each model. The parameters 
of Jhe models are the effective temperature T3 (in thousands of degrees), the 
logarithm of the surface gravity G, the mass fraction of hydrogen X, and the 
mass fraction of helium Y. The tabulated observables are as follows: UBV 
colors, uvby colors, bolometric corrections, Balmer discontinuity, helium 
discontinuity at X3422 ( 2 ^P-k), Balmer and Paschen slopes, and H 7 
equivalent widths and central intensities. 

The UBV colors were computed with the filter functions given by 
Matthews and Sandage (1963) and have been transformed, with their 
equations, to the standard Johnson and Morgan system. 

The Stromgren u-b and b-y colors and indices and have been 
computed with the filter functions of Graham (Matsushima, 1969). The 
normalization coefficients given by Matsushima (1969) were used to 
normalize them to the standard colors. The same normalization coefficients 
can be used since the “instrumental colors” of black bodies computed here 
are the same as Matsushima’s. Therefore, Matsushima’s coefficients put these 
models on the same system. 
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Table 1.— Opacity Sources 


Opacity 

Bound-bound hydrogen 
Neutral hydrogen. . . . 

Neutral helium 

Singly ionized helium . 
Negative ion of hydrogen 
Negative ion of helium . 


References 

Griem, 1960; Underhill, 1962 
'Mihalas, 1965 

Mihalas, 1965; Vardya, 1963 
Vardya, 1963 

Geltman, 1962; Fischel, 1963 

McDowell, Williamson, and 
Myerscough, 1966 


The bolometric correction is a color in the sense that it is the difference in 
magnitudes between two filters: 


BC = 2.5 log 





2.5 log 


r 




O 


After being normalized to the solar value (Allen, 1963), the resultant 
equation is 


BC = 12.92975 - 10 log Feff + 2.5 1ogjj^~ '^x’^X^^j 

where is the effective temperature of the model, is the flux per unit 
wavelength, and Vy^ is the visual filter function of Matthews and Sandage 
(1963). 

The Balmer discontinuity is computed as a magnitude difference at X3647 
between the flux per unit frequency shortward of X3647 and the flux F^ 
longward of X3647 extrapolated from X4586 and X4217 in a logarithmic plot 

F^(3647^) 

'^3647 ■ f^(3647-) 


where 


log Fj,(3647‘") = log F^(4217) + 


(3647 - 4217) 
(4585 - 4217) 


X 


log f 

F„(42I7) 


The wavelength X4217 is midway between and H^, and the wavelength 
X4586 is midway between and 
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The helium discontinuity is the magnitude difference between the Fj/s on 
either side of X3422: 


"^3422 


-2.5 log 


F^(3422^) 

F^(3422-) 


The Balmer and Paschen slopes are the rates of change of the magnitudes, 
with in the Balmer and Paschen continua: 

^ FJ3647*) / A 

4>d = -2.5 log AX = 0.0226/x 

^ ^ Fj,(3422) VAX/ 

and 


F„(5508) / 1 \ 

= -2.5 log (- 1 ), AX = 0.0458/x 

P Fj,(5050) VAX/ 

where AX is in microns. 

The profile of the has been computed with the ESW method (1967) 
for the Stark broadening, which was folded with a Voight function. In the 
solution of the transfer equation, electron scattering was included, and the 
total source function was allowed to vary across the line. The equivalent 
width was computed from 



The continuum flux F^ at each wavelength was found through the 
assumption that at ±70A the line flux F^ and the continuum flux are equal 
(i.e., the residual flux is zero) and then by using linear interpolation for 
F,(AX). 


TABLE OF THE MODELS 

Each pair of facing pages in Table 3 contains the data for one model. The 
left-hand page is comprised of the depth-dependent parameters. The column 
headings are- 

TAU ^5000 standard optical depth 

TE Temperature, °K 

PG Gas pressure, dynes/ cm^ 

PE Electron pressure, dynes/cm^ 

RHO Mass density, g/cm^ 

KAPZ Total opacity, absorption plus scattering at the 

^andard wavelength, cm^/g 
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DF(%) 

Percent error in the integrated flux 

DTE 

Temperature changes that would be applied to 
the next iteration, °K 

XH2 

Fraction .of hydrogen that is ionized 

XHE2 

Fraction of helium that is singly ionized 

XHE3 

Fraction of helium that is doubly ionized 


The right-hand page is the emergent flux (FLAM) in units of ergs/cm^/A as 
a function of wavelength (LAM) in angstroms. The values of X and Y are the 
mass fractions of hydrogen and helium. The effective temperature (r^^^ 
and logarithm of surface gravity (LOG G) are listed on each half of page. 
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Table 2 -Observables 
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